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ABSTRACT: We study at strong coupling the scaling function describing the large spin 
anomalous dimension of twist two operators in M = 4 super Yang-Mills theory. In the 
spirit of AdS/CFT duality it is possible to extract it from the string Bethe Ansatz equa- 
tions in the s[(2) sector of the AdS^ x S 5 superstring. To this aim, we present a detailed 
analysis of the Bethe equations by numerical and analytical methods. We recover several 
short string semiclassical results as a check. In the more difficult case of the long string 
limit providing the scaling function, we analyze the strong coupling version of the Eden- 
Staudacher equation, including the Arutyunov-Frolov-Staudacher phase. We prove that it 
admits a unique solution, at least in perturbation theory, leading to the correct prediction 
consistent with semiclassical string calculations. 

KEYWORDS: integrable quantum field theory, integrable spin chains (vertex models), 
quantum integrability (Bethe ansatz). 
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1. Introduction 

The n-gluon maximally helicity violating (MHV) amplitudes in planar J\[ = 4 SYM obey 
very remarkable iterative relations [Q, ||] suggesting solvability or even integrability of 
the maximally supersymmetric gauge theory. The main ingredient of the construction is 
the so-called scaling function defined in terms of the large spin anomalous dimension of 
leading twist operators in the gauge theory The scaling function can be obtained by 
considering operators in the s[(2) sector of the form 

Tr \v s Z L + permutations} . (1.1) 

The classical dimension is S + L, so L is the twist, with minimal value L = 2. The minimal 
anomalous dimension in this sector is predicted to scale at large spin S as 

A- S = f(g) log S + O(S ), (1.2) 

where the planar 't Hooft coupling is defined as usual by 

S 2 = ^, * = N 9 2 YM . (1.3) 
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The one and two loops explicit perturbative calculation of f(g) is described in [[| ^ ^] 
and H ^]. Based on the QCD calculation [|I]], the three-loop 0(g 6 ) calculation is per- 
formed in [ fiT| , |2p by exploiting the so-called trascendentality principle (KLOV). 

In principle, one would like to evaluate the scaling function, possibly at all loop order 
by Bethe Ansatz methods exploiting the conjectured integrability of SYM. This strategy 
has been started in [|l^]. In that paper, an integral equation providing f(g) is proposed 
by taking the large spin limit of the Bethe equations [14]. Its weak coupling expansion 
disagrees with the four loop contribution. The reason of this discrepancy is well under- 
stood. The Bethe Ansatz equations contain a scalar phase, the dressing factor, which is not 
constrained by the superconformal symmetry of the model. Its effects at weak-coupling 
show up precisely at the fourth loop order. 

A major advance was done by Beisert, Eden and Staudacher (BES) in [|l5p. In the 
spirit of AdS/ CFT duality, they considered the dressing factor at strong coupling. In that 
regime, it has been conjectured a complete asymptotic series for the dressing phase [jjj]]. 
This has been achieved by combining the tight constrains from integrability, explicit 1- 
loop cr-model calculations [ [P7] , [TJ, and crossing symmetry [p|[. By an impressive 

insight, BES proposed a weak-coupling all-order continuation of the dressing. Including it 
in the ES integral equation they obtained a new (BES) equation with a rather complicated 
kernel. The predicted analytic four-loop result agrees with the KLOV principle. Very re- 
markably, an explicit and independent perturbative 4-loop calculation of the scaling func- 
tion appeared in (j22[]. In the final stage, the 4-loop contribution is evaluated numerically 
with full agreement with the BES prediction. 

This important result is one of the main checks of AdS/ CFT duality. Indeed, a non 
trivial perturbative quantity is evaluated in the gauge theory by using in an essential way 
input data taken from the string side. 

As a further check, one would like to recover at strong coupling the asymptotic behav- 
ior of the scaling function, as predicted by the usual semiclassical expansions of spinning 
string solutions [[23[ f[4|, p>5| ]. Actually, the BES equation passes this check, partly numer- 
ically [ p6[ ] and partly by analytical means [|2^]. One could say that this is a check that 
nothing goes wrong if one performs the analytic continuation of the dressing phase from 
strong to weak coupling. 

From a different perspective, one would like to close this logical circle and check that 
the same result is obtained in the framework of the quantum string Bethe equations pro- 
posed originally in [17]. Indeed, it would be very nice to show that these equations repro- 
duce the scaling function in the suitable long string limit. Also, one expects to find some 
simplifications due to the fact that only the first terms in the strong coupling expression 
of the dressing must be dealt with. On the other side, the BES equation certainly requires 
all the weak-coupling terms if it has to be extrapolated at large coupling. 

In this paper, we pursue this approach. As a first step, we study numerically the 
quantum Bethe Ansatz equations in the 5 1(2) sector and check various results not directly 
related to the scaling function. Then, we work out the long string limit which is relevant to 
the calculation of f(g)- From our encouraging numerical results, we move to an analytical 
study of a new version of the BES equation suitable for the string coupling region. This 
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equation has been first derived by Eden and Staudacher in [ |13| ] as a minor result. Indeed, it 
has been left over because the main interest was focusing on matching the weak-coupling 
4-loop prediction. However, we believe that it is a quite comfortable tool if the purpose 
is that of reproducing the strong coupling behavior of the scaling function. We indeed 
prove that the solution described in p7[ ] is the unique solution of the strong coupling ES 
equation. 

The plan of the paper is the following. In Sec. @ we recall the Bethe Ansatz equations 
valid in the s[(2) sector of M = 4 SYM without and with dressing corrections. In Sec. @ we 
present various limits obtained in the semiclassical treatment of the AdS$ x S 5 superstring. 
We present our results for short and long string configurations. In Sec. @) we analyze the 
strong coupling ES equation building explicitly its solution and checking that it agrees 
with the result of [^]. We also investigate numerically the equation without making any 
strong coupling expansion to show that the equation is well-defined. Sec. (||) is devoted 
to a summary of the presented results. 

2. Gauge Bethe Ansatz predictions for the scaling function without dressing 

In the seminal paper JT3|], Eden and Staudacher (ES) proposed to study the scaling func- 
tion in the framework of the Bethe Ansatz for the s((2) sector of M = 4 SYM. The states to 
be considered in this rank-1 perturbatively closed sector take the general form 

Tr {{V s ' Z) ■■■ {V SL Z)}, Sl + --- + s L = S. (2.1) 

They are associated to the states of an integrable spin chain. The anomalous dimension A 
is related to the chain energies by 

A = L + S + E L , s {g) • (2.2) 

The all-loop conjectured Bethe Ansatz equations valid for E^^ up to wrapping terms are 
fully described in [^|[ Some explicit tests are can be found in |2^, ^(J. The Bethe 
Ansatz equations for the roots {v,k}i<k<s ar e 

9 2 

, . , r. „ , 1 - -F— 



^] = m: ! j - 4=x( Uk ± 1 -), (2.3) 



t + - t~ a 2 K V 2 



x k x j 



where we have defined the maps 



z(u) = |^l + yi-^j, u{x) = z + £. (2.4) 

The solutions of Eq. ( pT3| ) must obey the following constraint to properly represent single 
trace operators 

s x + 



EI ^ = !■ (2-5) 



fc=l X k 
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The quantum part of the anomalous dimension, i.e. the chain spectrum, is obtained from 

^si9) = 2^t{j,-j). (2.6) 

8=1 \ X k x k J 

Taking the large S limit of the Bethe Ansatz equations, ES obtain the following represen- 
tation of the scaling function 

f(g) = 8 5 2 - 64 5 4 r dta(t) (2.7) 

Jo 2gt 



where a is the solution of the integral equation 

t 



a{t) 



2gt 

with the (so-called main) kernel 



Wgt) _ Ag2 r dtlKm{2gt)2gtl)a[tl) 



(2.8) 



KmM= mMq^mm. (2 . 9) 

Notice that given a(t) we can simply write f(g) = 16 g 2 a(0) pl|]. These equations are 
independent on the twist which drops in the large S limit. This is important since the 
scaling function is expected to be universal [[32], ^] and therefore can be computed at 
large twist. 

Unfortunately, the perturbative expansion of f(g) disagrees at 4 loops with the ex- 
plicit calculation in the gauge theory. This is well known to be due to the missing contri- 
bution of the dressing phase. 

2.1 Input from string theory: Dressing corrections 

The effect of dressing is discussed in [|^] to which we defer the reader for general dis- 
cussions about its origin and necessity. The Bethe equations are corrected by a universal 
dressing phase according to 

1 

% = rr^^L — ^e 2 ^. (2.io) 

\x k J j&*t-*j 

x k x J 

The general perturbative expansion of the dressing phase is [0, [33p 1 

= E E E {9 2 Y +U+ " Pl? + \tl [qr(Pk)qr + i + 2v(Pj) - (* « j)] , (2-11) 

r>2 u>0 ii>v 



where the higher order charges q r (p) are [ p4| ] 



2sin(^p / Jl + 16 g 2 sin 2 f - 1 



r-1 



?r(p) = A 2 • P • ( 2 - 12 ) 

r — 1 l 4 g J sin § 1 



'This general formula holds unchanged in various deformations of the SYM theory see for ex. 
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(3) (2) 

The first non trivial constant is ^ 3 7^ 0- Indeed, P22, = consistently with the 3-loops 
agreement with explicit perturbation theory. 

The proposed coefficients for the all-order weak-coupling expansion of the dressing 
phase [fl6p are given in JIB] ] (see also [[3^] and [[3^]). They read 

« = 2( -ir- <^<^> („ _ 2 ;_+; + J ( 2 ;_ + ;)« 2 , + » 

and zero if/i — 7* — 1/ + 1 < 0. This proposal is completely equivalent to a precise modifi- 
cation of the kernel of the integral equation. It amounts to the replacement 

K m (t, t') -> K m (t, t') + 2 ff c (t, (2.14) 

where the dressing kernel can be written 

roo +11 

K c (t,t') = 4g 2 dt"K 1 {t,2gt")- [J ^—K (2gt",t'), (2.15) 
Jo e t - 1 

with 

KolMJ- <2 _ . ( 2 - i6 ) 

= ^7)2 • ( 2 - 17 ) 

The modified integral equation can be exploited to compute the perturbative expansion 
of /(<?). Now, there is agreement with the 4 loop explicit calculation. As we stressed in 
the Introduction, it is very remarkable that this weak coupling agreement is found with 
various inputs from string theory. In this sense, this is a powerful check of AdS / CFT 
duality. 

3. Strong coupling regime and the string Bethe equations 

As we explained, the BES equation is obtained by including in the ES equation an all- 
order weak-coupling expansion of the dressing phase. This expansion comes from a clever 
combination of string theory inputs and constraints from integrability. In our opinion, this 
is the essence of the integrability approach to AdS/ CFT duality. As a consistency check, 
one would like to recover from the BES equation the known semiclassical predictions 
valid in the AdS^ x S 5 superstring at large coupling. 

There are indeed several limits that can be computed. The semiclassical limit is eval- 
uated in terms of the BMN-like [f|(]] scaled variables which are kept fixed as A — > 00 

— — — (3D 

Vx' vT vT ( j 

where E, J,S are the semiclassical energy of a string rotating in S 5 with angular momen- 
tum J and spinning in AdS$ with spin S. The classical solution, and the first quantum 



corrections as well, are described in [23, 24, 23] 
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The simplest limits that can be considered are those describing short strings that do 
not probe AdS regions with large curvature. We call them 



short-GKP, < 1, < 1, : E = \/j 2 + 2SV\, (3.2) 

v A v A 

short-BMN, fixed, < 1, : E = J + S\ 1 + -^. (3.3) 

vA vA V ^ 



The scaling function is instead reproduced in the simplest long string limit which 
reads 

long string, -jj= < 1, -j= > 1, : £ = 5 + /(A) log 5. (3.4) 
In this limit, one can read the strong coupling behavior of the scaling function which is 



/string (5) =45- - log2 + --- . (3.5) 

7T 



A first attempt to solve numerically the integral equation is described in [^f]]. The 
BES equation is solved in a discrete series of Bessel modes and the result is a numerical 
profile of the scaling function shown in Fig. (1) of that paper. The lower curve is the 
scaling function taking into account the conjectured dressing and well reproduces the 
string prediction. 

Recently, the leading term of f(g) has been obtained analytically in p7|], including the 
analytic expression of the Bethe roots densities in the strong coupling limit. 



An intriguing feature of Fig. (1) in [ 26] is that the strong coupling behavior starts very 
early, namely at g ~ 1. The long string behavior should be visible at S/^/X 3> 1 which 
at g = 1 means S 3> 13. One can ask if it is possible to explore numerically the Bethe 
equations in the gauge theory up to g ~ 1 and with S = 0(100) to extract the scaling 
function. Actually, this is a hard task. At g ~ 1 the regime is not perturbative. The 
complete dressing should be resummed and it is not easy to do that, although some very 
interesting results have been presented in [JD-j]. 

An alternative hybrid approach would be that of taking the string Bethe equations 
with leading dressing. This should be enough to study the leading terms of f(g) at strong 
coupling. Of course, the problem is now that g must be large and then S must be unre- 
alistically large to deal with the numerical solution. However, not much is known about 
the properties of the strong coupling dressing expansion. It is divergent, but possibly 
asymptotic. Therefore, it would be difficult to estimate its accuracy at g ~ 1. 

In the next Section, we illustrate the detailed exploration of the above three limits. 



-6- 



3.1 The String Bethe equations and their numerical solution 

For the subsequent analysis it is convenient to pass to the p variables defining 

x + (u) 



p(u) 
x± (p) 



-ilog 



x~(u) ' 



+ 12 

e 2 
4sin§ 



1 -h ,/l + 16 5 2 sin 2 | j } 



u{p) = + sin 2 |. 



2 2 

The loop corrections to the energy are now 



fc=i 



^,s(5) = E(\/ 1 + 16 5 2sin2 y- 1 ). 



(3.6) 
(3.7) 
(3.8) 

(3.9) 



where {pk} are obtained by solving the Bethe equations with dressing Eq. Q2.10| ). At strong 
coupling, we use the leading dressing phase and write the string Bethe equations in loga- 
rithmic form as 



iLpk = E 

j^k 



( 



log 



\ 



x k x j 



J 1 

d 



+ 



(3.10) 



■ x j J 



+2i(u k - Uj)log 



X k X j 



X k X j 



\ 



x k x j 



> +2irie(p k ) 



(3.11) 



x k x j J 

where e(x) = x/\x\ and we have utilized the p — > — p symmetry of the solution for the 
ground state as well as its known mode numbers [13]. 

The numerical solution of the String Bethe equations is perfectly feasible. The tech- 
niques have already been illustrated in the two compact rank-1 subsectors su(2) and 



su(l|l) as discussed in gl], gj. 

First, we solve the equations at A = 0. This is the one-loop contribution. It is known 
exactly at L = 2. The Bethe roots are obtained, for an even spin S, as the roots of the 
resolvent polynomial [^, |J] 



n(tt - u s ) ~ 3 F 2 (-S, S + 1, - - i u; 1, 1; lj . (3.12) 

The 1-loop energy can also be computed in closed form with the result 

E 2 ,s = 8 9 2 (ip(S) - (1)) ~ 8 g 2 log S + > f(g) = 8 g 2 + ■ ■ ■ . (3.13) 

The roots are real and symmetrically distributed around zero. This features are also true 
for the ground state at arbitrary finite twist. 
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At L > 2, we use the solution at L = 2 as the starting point for the numerical root 
finder. Then, we increase g. At each step, we use a linear extrapolation of the previous 
solutions to improve the guess of the new solution. This procedure is quite stable and 
allows to explore a wide range of L, S, g values. Notice that changing the twist L is trivial 
since the complexity of the equations does not change. 

3.2 Short string in the GKP limit 

As a first numerical experiment we fix L and S and increase g up to large values where the 
equations are reliable. This is the short string limit where the AdS geometry is approx- 
imately flat. One expects to recover the Gubser-Klebanov-Polyakov law A ~ g 1 / 2 [pS]]. 
Solving the equations, we indeed verify that the Bethe momenta have the asymptotic 
scaling p k ~ g -1 / 2 . This is clearly illustrated in Fig. (||). 

Despite the non trivial distribution of the Bethe roots, it is straightforward to compute 
the anomalous dimension at large g. To do that, we first write the Bethe equations in the 
form 

iLp k = 2irie{ Pk ) + Y l Skj- (3.14) 

At large g we can write 

W = "?fe + — + ■■■• ( 3 - 15 ) 

g 1/2 9 

The Bethe momenta can be divided into a set P with a k > and a symmetric set obtained 
by flipping p — > —p. The expansion of the energy is 

E = E ( 2 9 1 ' 2 \0L k \ + 2 e(a k ) a' k - l) + ■ ■ ■ = (3.16) 
fc 

= 4 5 V2 ^ afc + 4 J>' fc -S. (3.17) 
fcep fcep 

From the Bethe equations we have 

keP kePjgp 

where we exploited Sij = —Sji. Each Bethe momentum pj with j & P can be written as 
—pji with j' G P. Hence we can write 

%L J2Pk=i*S+ Sfa,-Pj). (3.19) 

keP j,keP 

At large 3 

S{Pk, ~Pj) = -2t ^a fc ay + ^-(a4 a j + a fc a j)) H ( 3 -20) 

The leading order Bethe equations give 

iirS-2i oc k a J= — ► E a fc=(^) V2 - ( 3 - 21 ) 

fcjGP fcGP ^ ^ 
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The next-to-leading terms are 



%L a k = -2i ]T ( a 'k a j + a k a'j) = -4* a k a 'j, ( 3 - 22 ) 

keP j,keP keP jeP 

— E «'* = -f (3- 2 3) 



In conclusion, by plugging these results in the expression for A, Eq. (|L2|), we find 

A = v / 25A 1 / 4 + 0(A- 1 / 4 ). (3.24) 

The subleading term cancels the classical contribution leaving a pure A 1 / 4 behavior. This 
can be compared with the flat limit in the semiclassical approximation that reads for J/ vX 
5/VA « 1 



String = yJ 2 + 2SVX (3.25) 

with full agreement. 

3.3 Short string in the BMN limit 

If we keep S fixed and increase A with X/L 2 fixed, we can reach the BMN limit p0|]. 
This is numerically very easy because L enters trivially the equations. Fig. (Q) shows the 
convergence to the BMN limit when L is increased from 10 to 100 and S is fixed at S = 4. 
The various curves clearly approach a limiting one. This is very nice since it is an explicit 
show of how the BMN regimes sets up. Fig. (||) shows the limiting curves for S = 4, 6, 8 
at very large L = 10 4 . The three curves are perfectly fit by the expected law 

(3.26) 




3.4 Long string limit and the scaling function 

The previous pair of tests in the (easy) short string limits is a clear illustration that the 
numerical solution of the Bethe equations is reliable. 

The slow string limit is much more difficult. We begin with a plot of the energy at 
fixed twist L = 6 and increasing spin from 10 to 60. It is shown in Fig. (||). Each curve 
bends downward as g increases, since it ultimately must obey the g l l 2 law. However, at 
fixed g, when S increases the energy increases slowly eventually following the log S law. 
We attempted an extrapolation at S — > oo at each g. In Fig. (||) we show our estimate for 
the derivative of the scaling function, by fitting the data (S\ — S2) with Si < S < S2- We 
also show the analytical prediction 4. It seems to be roughly reproduced as soon as g > 1. 

The above "dirty" numerical procedure shows that it is reasonable to expect that the 
quantum string Bethe equations are able to capture the correct strong coupling behavior of 
the scaling function. However, the above extrapolation has a high degree of arbitrariness, 
especially concerning the fitting function employed to estimate the S — > 00 limit. Also, 
one would like to go to quite larger g requiring a huge number of Bethe roots, equal 
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indeed to the spin S. In practice, as it stands, the numerical investigation could hardly be 
significantly improved. 

For all these reasons, in the final part of this paper we explore an equation analogous 
to the BES equation, but derived for the string Bethe equation, at least with the leading 
order dressing phase. 



4. The strong coupling ES equation 

The inclusion of the AFS phase [ ]T7| ] in the ES equation 



o(t) 



e* - 1 



2gt 



4 5 2 / dt' K{2gt,2gt , )a{t') 



t-t' 



(4.1) 



(4.2) 



is straightforward and has been described in full details in [13]. The resulting equation 
reads 



a(t) 



t 



e* - 1 



Ji(2 5 t) , n J 2 {2gt) 



2gt 



2gt 



+ 2g ~ 4 ff 2 / dt 1 K(2gt,2gt')a(t l ) 



K(t,t') = K(t,t') + 2gK AFS (t,t'), 
where the AFS kernel is 



K t(J 2 (t)J (t')-J 2 (t')J (t)) 
K AF s(t,t) ~ ^T^7j2 ' 



, (4.3) 
(4.4) 

(4.5) 



It is important to remark that both the main and AFS kernels admit a simple expansion as 
series of products of two Bessel functions. For the main kernel, it is well known that 



(4.6) 



n>l 



For the AFS kernel we have 2 



J2{t)J {t')-J2{t')Mt) 8 „ 

1 2 _ U,)2 - l^ nJ 2n{t) J 2 n{t )■ 

* ' * ' n>l 

We change variables to put the equation in a somewhat simpler form and define 



(4.8) 



t 



2gJ 



2 This expansion can be proved by starting from the identity 



J 2 (*) Jo(i') - Mt 1 ) J (t) = — [t' Jx(i) J {t') - t Jo W Ji(t')] 



(4.9) 



(4.7) 



The result Eq. (4.8 ) follows immediately from the expansions reported in the appendices of [ |15| ] 



-10- 



The strong coupling ES equation for h is then 

Mr) = + 2 g - f dr> K(r, r<) ^ _ i Mr')- (4-10) 

This equation is expected to be reliable at strong coupling. It should reproduce the leading 
term in the scaling function f(g) = 8 g 2 h(0). With the conventions adopted in this part of 
the paper, this means h(0) = l/(2g) + subleading. 

4.1 Leading order at strong coupling 

We are interested in the large g limit with r fixed. The leading terms are obtained by 
expanding 

Mr) = ~*(r) + o(i), (4.11) 
- = 2g + 0(l). (4.12) 



e r/(2g) 

Taking the terms with the leading power of g we find that v(t) satisfies the remarkably 
simple equation 

H dr> K AFS (r, t 1 ) v(r') = (4.13) 
Jo 2.T 

Taking into account the expansion of the AFS kernel, the equation can be written in the 

equivalent form 



°° . J2n(r)v(r) 1 

dr - 2 = — S nil . (4.14) 



Following the approach of |26, 27] we expand the solution in a Neumann series of Bessel 
functions 

v(r) = Z^—> (4-15) 
fc>i r 

and now the equation reads 

T, H nkC k = ^5n,i, H nk = f°° J2n{X \ Mx) dx. (4.16) 

fe>l ib J ° x 

Now, the following question arises: Is this equation a constraint or does it determine a 
unique solution v(t)7 As a first step, we prove that the solution of [l27| ] is indeed a solution 
to the above equation. This solution reads 

+1 T(n + l/2) _ T(n + 3/2) 

C2n " r(n)r(i/2)' C2n+1 ~ r(n + i)r(i/2)- (417) 

Now, the detailed values of the matrix elements H n k are 

1 



k even : H ni k 



32 n(n + 1) (2n + 1) ! 



A; = 2n + 2, 



^ = 2n 

8n(2n + l)(2n-l)' ' (4.18) 

~ r r ~, ~ . Jc — 2tt< 2, 
32 (2n - l)n(n - 1) 

0, otherwise. 
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and 

16 (_i)«+p 
k = 1 + 2p, : H n , 1+2p = — _ — - 3)2)(4n2 _ (2p - 1)2 )( 4n 2 _ ( 2p _ 1)2) ( 4 - 19 ) 

The linear equations Eq. Q4.16D are then (£ stands for even Bessel) 

£ n = 0, n > 1, (4.20) 

where 

n-l n ' i; 32(2n - l)n(n - 1) 8n(2n + l)(2n - 1) 32n(n + l)(2n + 1) 

p>0 

These equations are indeed satisfied by the solution defined through Eq. ( |4.17| ). This can 
be checked by evaluating the infinite sum in closed form by the Sommerfeld-Watson trans- 
formation methods. For instance the five terms in Eq. ( |4.21| ) read at n = 1,2,3 

£ 1 = + — - — + 3L - — = 0, (4.22) 
48 256 768 16 

111 7 

6 2 = + = 0, (4.23) 

384 320 1024 15360 ' v ' 

115 311 , Mni . 

£ 3 = + + = 0. (4.24) 

1280 896 12288 430080 v ' 

However, this is not the unique solution of Eq. ( |4.13 ). As is quite usual, the straightfor- 
ward strong coupling limit of Bethe equations does not determine completely the solution 
which is fixed by the tower of subleading corrections. A similar difficulty is explained in 
details in [|7p. 

For instance, a second solution of Eq. ( [4.16D is 

c 2k = 2k{-l) k+ \ c 2fc+ i = 0. (4.25) 

Indeed, this produces the remarkably simple solution 

v{r) = £ 2M-l) fc+1 ^ = \ Mr). (4.26) 
fc>i r z 

Notice that, before summing the series, this second solution is precisely of the general 
class Eq. ( [4.15 ). As a check, we have indeed 

f°° , hn{r) Ji(t) 1 sin nir 1 

/ dT ^2 = ~o YT2 TT = T^ <5 ".i for n= 1,2,... (4.27) 

Jo 2t 2 8Trn 2 (n 2 -l) 16 

In practice, we still need the equal-weight condition on the even/ odd Bessel functions 
contained in the solution p7|]. 

To find a unique solution, we must examine the next orders in the strong coupling 
expansion. Indeed, the next orders provide both equations for the various subleading 
corrections to the solution v(t) and constraints on the previous contributions. This is due 
to the fact that the AFS kernel Kafs{t, r') is a function of r expressed as a Neumann series 
of purely even Bessel functions. The odd Bessel functions provide the above mentioned 
constraints as we now illustrate. 
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4.2 NLO order at strong coupling 

Let us work out the constraints from the subleading correction. If we take into account 
the next terms in the expansion and write 



h(r) = - V ^{t) + \v^{t) + --- 
9 r 

T 

2g 1 

3 2 



e r/{2g) _ i 

we find the following equation for 

/ dr' K afs (t,t')v^(t') = 
Jo 

= j~dr> ^K afs (t,t')t'-\k(t,t'^ „<«V) + ^. 



To simplify, we exploit 
1 



^K afs (t,t')t'--K(t,t , ) = 

2 1 

= — : Y] n J 2n {r) hnir 1 ) J~] n J n (r) J n [r') 

— f ^{2 n + 1) J 2 n+i(r) J 2n +i(T'). 



TT 



n>0 



Hence, the equation can be written 



HdT' K afs (t,t') «<V) = ±^l + J2(2n + l)p n 



where, explicitly, the constants @ n are given by 



G. 



n k 



L 



fc>i 

°° J 2n +i{x) J k {x) 

X 2 



dx. 



The relevant integrals are 



dx 



8n(2n + l)' 

8(n + l)(2n + 1) : 
0, 



k = n, 

k = n + 1, 
otherwise. 



and 



""nfe 



dx 



4(-l) 



n+fc 



TT (4 k 2 - (2 n + 1) 2 )(2 k + 2n- 1)(2 fc - 2 n - 3) 



(4.28) 
(4.29) 

(4.30) 
(4.31) 

(4.32) 
(4.33) 

(4.34) 



(4.35) 

(4.36) 
(4.37) 



(4.38) 



(4.39) 
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Equating to zero the coefficients of the odd Bessel functions J 2n +i(r) we obtain the con- 
straint on 

O n = 0, n > 0, (4.40) 

where 

O n = (1 - 5 n ) . ° 2n + + Y F nk c 2k -i ~ - S n o- (4.41) 

v ,; 8n(2n + l) 8(n + l)(2n + 1) ^ 4 n '° v 7 

These conditions are linear combinations of the previous equations. Indeed, one can check 
that 

S n = —(O n + O n _i). (4.42) 
4n 

So this constraint adds nothing new and, in particular, is satisfied by the Alday's solu- 
tion |]27||. Looking at the even Bessel functions, we obtain the homogeneous equation 

^■wf w <r = , (4 . 43) 







which admits the consistent solution i/ 1 ) = 0. 
4.3 NNLO order at strong coupling 

The next order in the \ j g expansion is what we need to fix uniquely v^°\ Expanding as 
before and writing now 

h{r) = -vW ( r ) + _1 V W ( T ) + A w (2) + . . . ( 4 .44) 
5 5 9 

2 

T 2 5 - - + ^- + ■ ■ ■ (4.45) 



e r/(2 5 )_ 1 * 2 243 
we find the following equation for t/ 2 ) 

4 / K"APs(r,r')v (2) (r')£fr' = (4.46) 
Jo 

= J™ [\r'K{r,r') - 1(t') 2 K aps (t,t')} v(°\r')dT' - „(°)(r)- 

The odd Bessel functions give a constraint independent on v^ 2 \ To compute it, we need 
the integrals (in our relevant range of values for n, k) 

r° J2 k+ i(x) J 2n (x) 2 (-i)- fc 

Jo x 7r(l + 2A:) 2 -4n 2 1 ; 

and 

[°° J 2k+ l{x) J 2 n+l{x) 1 , , 8 , 

y — z — dx= w^m n ' k - ( } 

One obtains immediately the crucial relation 

M-l rc2 I- 1 )"'' (449) 
2k + 1 7T 4- 2n (1 + 2&) 2 - 4n 2 ' K ' 

n>l v ' 
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This relation permits to write all odd coefficients C2k+i in terms of the even ones. Substi- 
tuting this relation in the truncated versions of the basic conditions 

£ n = 0, (4.50) 

one obtains a well-posed problem converging rapidly to the solution [ j27| ] without any a 
priori condition on the solution. For instance, by truncating the problem dropping all C2 n 
with n > N, we find the table Tab. (fl|) of values for C2. A simple polynomial extrapolation 
to N — » oo provides the correct limit c\ = C2 = 0.5000(1). 

Hence, the strong coupling expansion is well defined and the leading solution is 
unique. Of course, it is the one described in [^]. 



N 


10 


20 


30 


40 


50 


60 


70 


80 


C2 


0.5665 


0.5386 


0.5273 


0.5211 


0.5173 


0.5146 


0.5127 


0.5112 



Table 1: Coefficient C2 from the truncated full rank linear problem. 



4.4 Numerical integration of the strong coupling ES equation 

To summarize, we have shown that the strong coupling ES equation is consistent with 
the results of [27]. In that paper, it was crucial to fix the relative weights of the even/odd 
Bessel functions appearing in the general solution. These weights were shown to be more 
than an Ansatz. They are encoded in the full equation before expanding at strong cou- 
pling. Alternatively, they can be derived by analyzing the next-to-leading and next 2 -to- 
leading corrections. 

As a final calculation and check, we provide the results from a numerical investigation 
without any strong coupling expansion to see how the correct strong coupling solution 
arises. This can be done along the lines illustrated in [26, ^7] ]. We start again from the 
Neumann expansion 

fc(T) = 5>— , (4.51) 
fc>i r 

and arrive at the infinite dimensional linear problem 

x — ^ m f 0, nodd 

c n = 6 nA + 2g5 ni2 -2n}2Klc k -8g \ _ ^ ^ ^ _ _ (4.52) 

fe>i 

where 




neven 



z M {g) = / ^p^ _ (4 . 53) 

As in [26, [2^] we can truncate this equation setting c k = for k > N. The solution should 
be reliable for g < N. We considered 2 < g < 30 and N = 50. The best fit to the leading 
term in f(g) gives 

f(g) = eg + subleading, (4.54) 
c = 4.000(1). (4.55) 

This confirms that the strong coupling ES equation has a unique solution with the correct 
large g limit. Of course, as g — > oo it is nothing but Eq. ( |4.17| ). 
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5. Conclusions 



In this paper, we have considered several properties of the quantum string Bethe equa- 
tions in the st(2) sector with the leading strong coupling dressing, i.e. the AFS phase. We 
have performed a numerical investigation of the equations showing that their analysis 
is quite feasible. As an interesting result, we have repeated the calculation of the GKP 
limit of the anomalous dimensions as for the highest excited states in the compact rank-1 
subsectors su(2) and su(l|l). Also, we have been able to observe the setting of the BMN 
scaling regime by reproducing the plane wave energy formula at fixed spin and large 
twist. In the case of the long string regime, we have been able to provide numerical ev- 
idences for a scaling function exhibiting an early strong coupling behavior as expected 
from the numerical solution of the BES equation. 

Motivated by these results, we have analyzed analytically and perturbatively at strong 
coupling an almost trivially modified version of the BES equation with the very simple 
strong coupling dressing [O]. In particular, we have proved that this equation admits, as 
it should, a unique solution for the asymptotic Bethe root (Fourier transformed) density 
in full agreement with existing results. 

While this work was under completion, the paper [ fi5l ] presented an analysis partially 
overlapping with our results. That paper derives an integral equation for the Bethe root 
density taking into account the dressing at strong coupling and is based on a novel integral 
representation of the dressing kernel. We hope that the two alternative approaches will 
turn out to be useful in computing the one-loop string correction to the large g scaling 
function. Indeed, this interesting contribution has been checked numerically in p6| ] but it 
still evades an analytical confirmation. 

Hopefully, these various efforts might give insight on the general structure of the 
dressing phase as well as on the role of the asymptotic Bethe equations in an exact de- 
scription of the planar spectrum [|46|]. Significative studies about finite size effects p7| ] 
and corrections that arise in a finite volume to the magnon dispersion relation at strong 



coupling Q48I] , see also [|49p, as well as the recent observation [ pOQ that the dressing phase 
could originate from the elimination of "novel" Bethe roots, strongly demand a deeper 
understanding . 
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Figure 1: Numerical solution in the short-GKP limit, i.e. at L, S fixed and increasing g. We show 
the p -1 / 2 behaviour of the Bethe momenta, the computed energy with its g 1 / 2 profile, and the 
(vanishing) difference with analytic strong coupling expansion Eq. (3.24). 
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S = 4 




Figure 2: Short-BMN limit. Convergence to the BMN limit at L, A — » oo with A/L 2 fixed. 
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Figure 3: Short-BMN limit. Check of the semiclassical string prediction String = L + S^/l + \/L 2 . 
We plot the quantum part E L<S computed by the Bethe Ansatz. 
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Figure 5: Long string limit. Best fit results for the asymptotic derivative of the scaling function. 
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